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The LHC has already provided many relevant measurements for the determination of parton 
distribution functions (PDFs). Measurements of the W and Z lepton distributions are of 
interest for flavor separation and in particular for the determination of the relatively poorly 
constrained strange quark distribution. In this contribution we shall discuss the computational 
developments that allow for the efficient inclusion of LHC data into the NNPDF framework 
consistently at NLO for all observables, and we study the constraints of the LHC W and Z 
data on the strangeness content of the proton. 



There have been a number of experimental measurements of direct relevance to PDF de- 
termination performed by LHC collaborations. Me asurem ents of inclusive jet and dijet cross 
sectiomPSE] anc j electroweak vector boson productioA^^ provide information on PDFs in pre- 
viously unexplored kinematical regions. While the importance of including LHC data in future 
determinations is clear, of particular interest is the potential impact of these data sets upon 
collider only fits, where existing determinations tend to be poorly constrained. This necessi- 
tates the inclusion into PDF determinations low energy data with potential contamination from 
nuclear corrections or higher twist effects. 

However, including collider measurements into a PDF fit on a large scale requires substantial 
computational resources. NNPDF parton sets are fitted by genetic algorithm minimisation over 
a large number of generations^ therefore the NNPDF methodology requires a fast method of 
computing collider observables. Although LHC data was previously included in NNPDF2.iPby 
a reweighting method, the constraining power of the LHC dataset makes adding a large quantity 
of data in this manner impractical. We shall here describe a fast convolution method that has 
been developed to enable the inclusion of new hadronic data in NNPDF fits, before going on 
to discuss some preliminary results from fits including LHC data (NNPDF2.3 preliminary). 
There has been particular interest in the usefulness of the recent ATLAS measurements of W/Z 
productiorP in providing information on the strange content of the protonP We shall discuss 
here some preliminary results on the proton strangeness fraction using the updated fit. 

A number of tools are available for the computation of hadronic observables that allow for 
a straightforward variation of the input PDF a posteriori, a prerequisite for utility in parton 
fitting. In particular, the FastNL(P2 an( j APPLGriJP projects provide software which is well 
suited for use in fitting. The principle of these projects is to store the required perturbative 
coefficients for a process as weights upon an interpolating grid in x and Q space. The convolution 
required to calculate the observable is then reduced to a simple product, the PDF in the product 
may be straightforwardly varied along with the chosen value of as- For example, to compute a 



hadronic cross section in the APPLGrid framework, the following calculation is performed, 
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where the indices a, /3 run over points in the x-space grid, r runs over points in Q 2 , p denotes the 
perturbative order of the contribution, and I denotes the specific parton level subprocess. The 
W table contains the values of the Monte Carlo weights for a particular subprocess point, and 
the F® are the incoming subprocess parton densities constructed as a combination of PDFs as 
appropriate for the process in consideration. This method of computing observables is fast, but 
substantial speed improvements can be gained by combining PDF evolution with this procedure. 
For a set of flavour basis PDFs /, we write a general subprocess density as, 

13 

F®(x a ,x p ,Cfi) =J2C% (Mx^QDMx^Ql)) . (2) 

Where i,j denote the PDF flavour, and the cf^ are coefficients specifying how the subprocess 
density I is to be built. The evolution of the initial state PDF to the required scale Q 2 can 
be performed in an analogous fashion to the convolution in Eqn[T]by evaluating the matrix of 
DGLAP evolution kernels upon an interpolation grid as per the FastKernel method^ Obtaining 
the evolved PDF is reduced once again to a product, 
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Where the N are PDFs in a suitable evolution basis that diagonalises the matrix of DGLAP evo- 
lution kernels. The matrix EJ k holds the values of the DGLAP evolution kernel Tj^ (x,Qq,Q 2 ) 

convoluted with the interpolating basis functions as in reference^, and the matrix R is the 
rotation matrix from the evolution to the flavour basis. Here we adopt the notation N° for the 
-^pdf light evolution basis PDFs parameterised at the initial fitting scale Qq. It is now simple 
to construct the subprocess density using these matrices, 
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With the PDF evolution factorized, the computation in Eqn[T]is now reduced to a much simpler 
form particularly suited to a fitting application, 
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is the weight matrix containing all the values that may be precomputed and stored prior to a 
PDF fit. The calculation of a hadronic observable is then simply a matter of a sum of products 
over a grid in x-space, and the now reduced flavour basis of iVpdf light PDFs. Through this 
method we are able to reproduce the results of the original APPLGrid/FastNLO calculation at 
the same level of precision and with a substantial improvement in speed. 



Using this technique, we can now present results on the strangeness fraction with recent 
preliminary NNPDF fits including LHC data. R s (x,Q 2 ) = (s + s)/(u + d) has been deter- 
mined from a number of NNPDF fits at NLO to different datasets. Firstly a fit (here denoted 
NNPDF2.3 prelim) to the full NNPDF2.1 dataset with the addition of ATLAS 35 pb" 1 inclusive 
jet measurements^ ATLAS 35 pb _1 W and Z rapidity distributions^ and CMS 840 pb _1 W 
electron asymmetry dataff^ Secondly, a fit exclusively to the NNPDF2.3 collider data subset 
(NNPDF2.3 Collider), and finally, a fit to the HERA-I combined datasePS and ATLAS W/Z 
measurement only (NNPDF2.3 HERA+ATLASWZ). For comparison, the value of R s (x,Q 2 ) 
determined by the NNPDF2.t^ fit is also provided. 

A study by the ATLAS collaboration on the strange content of the protonP based upon 
fits to the same dataset as NNPDF2.3 HERA+ATLASWZ suggest that the ratio of strange 
to non-strange PDFs may be underestimated by previous determinations from global fits. In 
Figure Q] we examine the ratio of strange to non strange PDFs for the NLO fits NNPDF2.3 
prelim, NNPDF2.3 HERA+ATLASWZ, and NNPDF2.1. From this figure it is clear that the 
recent ATLAS W/Z measurements provide a valuable constraint, however at medium to large- 
x the HERA and ATLAS data alone is insufficient to provide a precise determination of the 
strangeness. 
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Figure 1: Comparison of the proton strangeness fraction determined from fits to various datasets. 

In Table 1 we compare the values obtained by the different sets at specific values of (x,Q 2 ) 
and see a similar pattern. Fits to reduced datasets, such as the collider only and HERA + 
ATLAS W/Z fits suggest a higher value of R s , however the data provides little constraint and 
therefore the uncertainties are substantially larger than in the determinations provided by the 
global fits. The values all broadly agree within the large uncertainties of the HERA + ATLAS 
W/Z fit as shown in the comparison in Figure El We can therefore conclude that the collider 
data alone is not yet sufficiently constraining to provide a precise determination of the proton 
strangeness fraction, and that the uncertainty on R s in the ATLAS determinatioiP has been 
underestimated. 

Acknowledgments 

The author would like to thank the organisers of the Rencontres de Moriond 2012 for the 
opportunity to present this work. This contribution made use of resources provided by the 
Edinburgh Compute and Data Facility (ECDF) ( |ht t p : / / www .ecdf.ed.ac.uk/[ ). The ECDF is 
partially supported by the eDIKT initiative (http://www.edikt.org.uk). 



PDF Set 


i? s (0.013, Mj?) 


R s (0.023, 1.9GeV a ) 


NNPDF2.1 
NNPDF2.3 preliminary 


0.61 ±0.09 
0.68 ± 0.06 


0.24 ±0.09 
0.36 ±0.10 


NNPDF2.3 HERA+ATLAS WZ only 
NNPDF2.3 Collider only 


1.00 ± 0.33 
1.00 ±0.28 


1.40 ±2.20 
0.95 ± 0.60 



Table 1: Table of R s values determined from several PDF fits and at two choices of (x,Q 2 ) 

Strangeness Fraction R (Q 2 = M 2 x = 0.013) Strangeness Fraction R (Q 2 = 1 .9 GeV 2 , x = 0.023) 
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Figure 2: Value of R s determined with different PDF sets at two choices of (x, Q 2 ). 
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